The Solar hep Process in Effective Field Theory 
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Using effective field theory, we calculate the S'-factor for the hep process in a totally parameter- 
free formulation. The transition operators are organized according to chiral counting, and their 
matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated- 
hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy- 
baryon chiral perturbation theory are considered. Fixing the only parameter in the theory by fitting 
the tritium /3-decay rate, we predict the hep S'-factor with accuracy better than ~ 20 %. 
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The hep process 

p + 3 Hc 



l He- 
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figures importantly in astrophysics and particle physics; 
it has much bearing upon issues of great current inter- 
est such as the solar neutrino problem, non-standard 
physics in the neutrino sector, etc. The hep reaction 
produces the highest-energy solar neutrinos even though 
their flux is much smaller than that of the 8 B neutrinos. 
At some level, therefore, there can be a significant distor- 
tion of the higher end of the 8 B neutrino spectrum due 
to the hep neutrinos. This change can influence the in- 
terpretation of the results of a recent Super-Kamiokandc 
experiment that have generated many controversies re- 
lated to the solar neutrino problem and neutrino oscil- 
lations |l|,|[. To address these issues, a reliable estimate 
of the hep cross section is indispensable. Its accurate 
evaluation, however, presents a challenge for nuclear and 
hadron physics [|| . The degree of this difficulty may be 
appreciated by noting that theoretical estimates of the 
hep astrophysical 5-factor have varied by orders of mag- 
nitude in the literature. The reason for this great vari- 
ation is multifold. Although the primary hep amplitude 
is formally of the Gamow- Teller (GT) type (A J = 1, 
same parity), the single-particle GT matrix element for 
the hep process is strongly suppressed due to the sym- 
metries of the initial and final state wave functions. Fur- 
thermore, the main two-body corrections to the "lead- 
ing" one-body GT term generically come with opposite 
sign causing a large cancellation. Therefore, it is neces- 
sary to calculate these corrections with great accuracy, 
which is a highly non-trivial task. The nature of the 
specific challenge involved here can be further elucidated 
in terms of the chiral filter picture. If meson-exchange 
current (MEC) contributions to a specific transition am- 



plitude are dominated by the one-soft-pion exchange di- 
agram, then one can take advantage of the fact that the 
soft-pion-exchange piece is uniquely dictated by chiral 
symmetry and that there is a mechanism (called the chi- 
ral filter mechanism) that suppresses higher chiral-order 
terms [0J|] . We refer to a transition amplitude to which 
the chiral filter mechanism is applicable (not applica- 
ble) as a chiral-protected (chiral-unprotected) case. It 
is known that the space component of the vector cur- 
rent and axial-charge are chiral-protected, whereas the 
GT transition is chiral-unprotected. This means that, 
in contrast to isovector Ml and axial charge transitions, 
both of which are chiral-protected and hence easily cal- 
culated Jg,0|, the primary hep amplitude, i.e. the GT 
transition, needs to be calculated without invoking the 
chiral-filter mechanism and hence is more subtle. 

In a highly successful method, which we call here the 
standard nuclear physics approach (SNPA) g, one uses 
the potential picture supplemented with MEC contribu- 
tions, the structures of which are obtained from one- 
boson- exchange. Many applications of SNPA to elec- 
troweak transitions in few-body systems are well doc- 
umented. A detailed SNPA calculation for the hep cross 
section was recently carried out by Marcucci et al. [||; 
this study has re-confirmed the substantial cancellation 
between the one-body and two-body terms for the GT 
transition p0| , pd| . It is therefore of primary impor- 
tance to have a good theoretical control of short-distance 
physics. A "first-principle" approach based on effective 
field theory (EFT) will hopefully provide an insight into 
this issue. A possible approach that is formally consis- 
tent with systematic power counting is the pionlcss EFT 
(for a recent review, see [p2|). This approach, however, 
cannot be readily extended to systems with A > 3. Apart 
from the basic problem of organizing chiral expansion for 



complex nuclei from "first-principles" , a plethora of pa- 
rameters involved would present a major obstacle. This 
difficulty is expected to be particularly pronounced for 
the hep reaction. 

In this paper we calculate the hep astrophysical S- 
factor, adopting a variant of EFT, that has been proven 
to be extremely successful in estimating the S'-factor for 
solar pp fusion p3|. The same method has also scored 
great success in predicting polarization observables in 
thermal np capture |l4|| . Our starting point is the ob- 
servation that, to high accuracy, the leading-order single- 
particle matrix elements in SNPA and EFT are identical, 
and that they can be reliably estimated with the use of 
realistic SNPA wave functions for the initial and final 
nuclear states. Next, we note that in EFT the operators 
representing two-body corrections to the leading-order 
one-body term can be controlled by systematic chiral ex- 
pansion in heavy-baryon chiral perturbation theory ]15[ . 
Then, since the ratio of a two-body matrix element to 
the leading-order one-body matrix element can be evalu- 
ated with sufficient accuracy with the use of the realistic 
SNPA wave functions, we are in a position to obtain a 
reliable estimate of the total (one-body + two-body) con- 
tribution. This method, which exploits the powers of both 
SNPA and EFT, may be referred to as a "more effective 
EFT (MEEFT)". Here we shall work out MEEFT up to 
next-to-next-to-next-to-leading order (N 3 LO), adopting 
a cut-off regularization. Our method, which to this or- 
der preserves chiral symmetry and renormalization group 
invariance, enables us to predict the hep S'-factor in a to- 
tally parameter-free manner. 

In the present scheme it is sufficient to focus on "irre- 
ducible graphs" in Weinberg's classification [[16| . Graphs 
are classified by the chiral power index v given by v — 
2(j4 — C) + 2L + J2i v ii where A is the number of nu- 
cleons involved in the process (A=4 in our case), C 
the number of disconnected parts, and L the number of 
loops. The chiral index, i/$, of the i-th vertex is given 
by Vi = di + e, + rij/2 — 2, where d,, e, and n% are re- 
spectively the numbers of derivatives, external fields and 
nucleon lines belonging to the vertex. The Feynman di- 
agrams with a chiral index v are suppressed by {Q / K x ) u 
compared with the leading-order one-body GT operator, 
with Q standing for the typical three-momentum scale 
and/or the pion mass, and A x ~ tojv ~ 47r/ w for the chi- 
ral scale. The physical amplitude is then expanded with 
respect to v. 

In this paper we shall limit ourselves to N 3 LO, al- 
though it is possible to go to N 4 LO without introducing 
new parameters. We write the current as 



J"(g) = V»(g) + A»{q) = f dx c - zqx J^(x), 



(2) 



where q is the momentum carried by the lepton pair. 
The calculation is considerably simplified by the facts 
that we are dealing with small momentum transfers and 



that, as a consequence, both q = \q\ and the time com- 
ponent of the nucleon momentum are 0(Q 2 /A x ) rather 
than 0(Q) as naive counting would suggest. For the 
vector and axial one-body current and charge operators, 
and also for the vector-current and axial-charge (chiral- 
protected) two-body operators, we can simply carry over 
the expressions given in Ref. H), hereafter referred to 
as MSVKRB. The vector charge two-body operator does 
not appear up to the order considered, while the EFT 
axial-current two-body operator is given, in momentum 
space, in Ref. 13 1. These two-body currents are valid 
only up to a cutoff A. This implies that, when we go to 
coordinate space, the currents must be regulated. This 
is a key point in our approach. We introduce a Gaus- 
sian cutoff regulator in performing Fourier transforma- 
tion. The resulting r-space expressions of the currents in 
the center-of-mass (cm.) frame are 



V 12 (r) = 

A? 2 (r) = 

A 12 (r) = 



9A m l 



12 f2 T x V [""1 ' °"2 2/oaM + Sl2 S/SaM] 



A 



q x 



OxJ/SaM + ( T x - -e» x ) yUr) 



-q ■ r cr ■ r 



9A_ - 

9Aml 
'2m N f^ 



ViaW, 



+ 



e 3 (T + + r_) 



c 4 + - I T x 



(J A 



2m N fl 



-r x {p 1 cr 2 ■r+p 2 cr 1 r 



2/oaW 

VIST) 



+5 A (r)d R O x , 



where r = r\ — r 2 , Si 2 is the tensor operator, and 



(3) 



Mr), yZ A (r)} 



dk 

(27T) 5 



e -fc7A 2 e »fc-r 



1, 



1 



k 2 



2/iaM = -r[y% A (r)}' and y% A (r) = (r/ml)[[y^ A (r)}'/r}' 
(the ' symbols denote derivatives), 0\ 



T cr ©> 



Or. 



T Q cr 0l T = rr-e> Q -±e> ,© = ±, x,r ee (ti0t 2 ) = 
(n © t 2 ) x - i(n t 2 ) v and (T ee (<n <r 2 ). Finally, the 
parameter d R is given by 



d H = d x + 2d 2 



1 



1 



3 C3 + 3 C4 + 6 



(4) 



The dimensionless parameters, c's and eTs, are defined 
in Ref. [0. The numerical values of C3 and £4 have 
been determined previously [[L8| : C3 = —3.66 ± 0.08 and 
£4 = 2.11±0.08. The derivative operators, pj ee h(p'i+Pi) 
(i = 1, 2) in Eq.(||), should be understood to act only on 
the wave functions. 



The explicit degrees of freedom in our scheme are the 
nucleon and the pion, with all other degrees of freedom 
(p- and w-mesons, A(1232), etc.) integrated out. There- 
fore, a reasonable range of the cutoff A would be some- 
where around 500-800 MeV. We shall consider here three 
exemplary values, A = 500, 600, 800 MeV. 

A crucial point is that Avi in Eq.(0) contains only 
one unknown parameter, d R , that needs to be fixed us- 
ing an empirical input. The tritium /3-decay rate, Tp, 
can be used for this purpose, since it is dominated by 
the unsuppressed GT term [ [19] . Thus, for each value 
of A, we adjust d R to reproduce the experimental value 
of r^ and, using the value of d R so determined, we 
evaluate the hep amplitude. To carry out this pro- 
gram, we must calculate both the tritium /3-decay am- 
plitude and the hep amplitude for the above-described 
currents, using realistic wave functions. We adopt here 
the method employed by MSVKRB. In particular, we 
use the correlated-hyperspherical-harmonics wave func- 
tions [g0|, 21 1 obtained with the Argonne Vis (AVI 8) two- 
nucleon [22] and Urbana-IX three- nucleon [£3| interac- 
tions. 

In the notation of MSVKRB, the GT-amplitudes are 
given in terms of the reduced-matrix elements (RMEs) 
Li(q;A) and E\(q;A). Since these RMEs are related to 
each other as Ei(q;A) ~ \/2Zi(g; A), with the exact 
equality holding at q=0, we consider here only one of 
them, Li(q; A). For the three exemplary values of A, Ta- 
ble I gives the corresponding values of d R , as determined 
from r^ |n|], and L\{q; A) at g=19.2 MeV and zero cm. 
energy. We see from the table that the variation of the 
two-body GT amplitude (row labelled "2B-total" ) is only 
~10 % for the range of A under study. The A-dependence 
in the total GT amplitude is made more pronounced by a 
strong cancellation between the one-body and two-body 
terms, but this amplified A-dependence still lies within 
acceptable levels. 



A (MeV) 


500 


600 


800 


d R 


1.00 ±0.07 


1.78 ±0.08 


3.90 ±0.10 


Life A) 


-0.032 


-0.029 


-0.022 


IB 


-0.081 


-0.081 


-0.081 


2B (without d R ) 
2B (ex d R ) 


0.093 
-0.044 


0.122 
-0.070 


0.166 
-0.107 


2B-total 


0.049 


0.052 


0.059 



TABLE I. Values of d R and Li(q; A) (in fm 3/2 ) calculated 
as functions of the cutoff A. The individual contributions 
from the one-body (IB) and two-body (2B) operators are also 
listed. 

Table [jl] shows the contribution to the S-factor, at 
zero cm. energy, from each initial channel. For com- 
parison we have also listed the MSVKRB results for the 
AV18/ULX interaction. It is noteworthy that for all the 
channels other than 3 Si, the A-dependence is very small 



A (MeV) 


500 


600 


800 


MSVKRB 


L S 


0.02 


0.02 


0.02 


0.02 


3 Si 


7.00 


6.37 


4.30 


6.38 


3 Po 


0.67 


0.66 


0.66 


0.82 


'Pi 


0.85 


0.88 


0.91 


1.00 


3 A 


0.34 


0.34 


0.34 


0.30 


3 P 2 


1.06 


1.06 


1.06 


0.97 


Total 


9.95 


9.37 


7.32 


9.64 



TABLE II. Contributions to the S'-factor (in 10~ 20 keV-b) 
from individual initial channels calculated as functions of A. 
The last column gives the results obtained in MSVKRB. 



(< 2 %). While the GT terms are dominant, the con- 
tribution of the axial-charge term in the 3 Si channel is 
sizable even though it is kinematically suppressed by the 
factor q. It is therefore reassuring that the chiral-filtcr 
mechanism allows a reliable evaluation of this amplitude. 
Summarizing the results given in Table O, we arrive at 
a quantitative prediction for the hep S'-factor: 



5*= (8.6 ±1.3) x 10- 20 kcV-b, 



(5) 



where the "error" spans the range of the A-dependcncc 
for A=500-800 MeV. This result should be compared to 
that obtained by MSVKRB ^S = 9.64 x lO" 20 keV-b. 
Note that the earlier studies |l0|]ll|] were based on less 
accurate variational wave functions than used here and 
in MSVKRB, and did not include P-wave capture contri- 
butions, which account for ~ 40 % of the total S-factor. 
To decrease the uncertainty in Eq. (g) , we need to reduce 
the A-dependence in the two-body GT term. According 
to a general tenet of EFT, the A-dependence should di- 
minish as we include higher order terms. A preliminary 
study Q indicates that it is indeed possible to reduce 
the A-dependence significantly by including N 4 LO cor- 
rections. 

In order to better understand how the present scheme 
works, it is helpful to compare the hep reaction with 
the radiative np-capture. The polarization observables in 
ft + p — ► d ± 7 are known to be sensitive to the isoscalar 
Ml matrix element, MIS', and this amplitude has been 
extensively studied in EFT J14U23]. The similar features 
of the hep GT amplitude and the MIS matrix element 
are: (i) the leading one-body contribution is suppressed 
by the symmetries of the wave functions; (ii) there is no 
soft-pion exchange contribution; (iii) nonetheless, short- 
range physics can be reliably subsumed into a single con- 
tact term. In the ftp case the strength of this term can be 
determined from the deuteron magnetic moment (for a 
given value of the cutoff A) . The calculation in Ref . JlJ] 
demonstrates that the A-dependence in the contact term 
and that of the remaining terms compensate each other 
so that the total MIS is stable against changes in A. This 
suggests that, if we go to higher orders, the coefficient of 
the contact term in question will be modified, with part 



of its strength shifted to higher order terms; however, 
the total physical amplitude will remain essentially un- 
changed. These features are quite similar to what we 
have found here for the hep GT amplitude. 

Evaluating the matrix element of the leading-order 
one-body operator in EFT with the use of realistic nu- 
clear wave functions is analogous to fixing parameters 
in an EFT Lagrangian (at a given order) using empir- 
ical inputs. The realistic wave functions in SNPA can 
be regarded as a theoretical input that fits certain sets 
of observables [£6|. In the present MEEFT scheme, we 
take the view that the same realistic wave functions also 
provide a framework for reliably calculating corrections 
to the leading-order one-body matrix element. While 
from a formal point of view the approach adopted here 
is, in certain cases, not in strict accordance with the sys- 
tematic power-counting scheme of EFT proper, never- 
theless the severity of this potential short-coming may 
depend on individual cases (see discussion in |2j). For 
the hep amplitude under consideration here, the mod- 
erate A-dependence exhibited by the numerical results 
suggests that the lack of rigorous power-counting cannot 
be too significative. Indeed, this type of "resilience" may 
also explain why the SNPA calculation in Ref. M gives a 
result very similar to the present one. It is true that the 
MSVKRB two-body terms are not entirely in conformity 
with the chiral counting scheme we are using here; some 
terms corresponding to chiral orders higher than N 3 LO 
are included, while some others which are N 3 LO in EFT 
are missing. This formal problem, however, seems to be 
largely overcome by the fact that also in MSVKRB a pa- 
rameter (the axial ttNA coupling strength) is adjusted 
to reproduce Tp. 

The latest analysis of the Super-Kamiokande data |2q | 
gives an upper limit of the solar hep neutrino flux, 
<&(hep) SK < 40 x 10 3 cm _2 s _1 . The standard solar model 
p9| using the hep S'-factor of Marcucci et al. |9[ predicts 
<£>(hep) SSM = 9.4 x 10 3 cm^s" 1 . The use of the central 
value of our estimate, Eq.(||), of the hep S'-factor would 
slightly lower $>(hep) but with the upper limit com- 
patible with $(/iep) SSM in Ref. |2J|. A concrete estimate 
of the theoretical uncertainty in Eq. (H) is expected to be 
useful for further discussion of the solar hep problem. 
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